ESD-TR-67-393 

ESTI  FILE  COPY 


CSD  RECORD  COPY 

RETURN  TO 

ESD-TR-67-393  SCIENTIFIC  8  TFCHNIf.  ,  INf-vi.RWATION  lMS'iE'N-66 

ESD  ACCESSION  LIST  81"‘D™e  >*n 

ESTI  Call  No*  AL  578  5.’  _ 

Copy  No  _ f_  of  ^  ,  oya. 

ERROR  PERFORMANCE  OF  DIFFERENTIAL  PHASE  SHIFT 
TRANSMISSION  OVER  A  TELEPHONE  LINE 

I 

AUGUST  1967 
J.  J.  Bussgang 
M.  Leiter 


Prepared  for 

DEPUTY  FOR  SURVEILLANCE  AND  CONTROL  SYSTEMS 
AEROSPACE  INSTRUMENTATION  PROGRAM  OFFICE 

ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
L.  G.  Hanscom  Field,  Bedford,  Massachusetts 


i 

I 

I 

J 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  un¬ 
limited. 


Project  705B 

Prepared  by 

TIIE  MITRE  CORPORATION 
Bedford,  Massachusetts 

Contract  AF  19(62S)-5  1 65 


When  US  Government  drawings,  specifications,  or 
other  data  are  used  for  any  purpose  other  than  a 
definitely  related  government  procurement  operation, 
the  government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications, 
or  other  data  is  not  to  be  regarded  by  implication 
or  otherwise,  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related 
thereto. 


Do  not  return  this  copy.  Retain  or  destroy. 


ESD-TR-67-393 


MTP-66 


ERROR  PERFORMANCE  OF  DIFFERENTIAL  PHASE  SHIFT 
TRANSMISSION  OVER  A  TELEPHONE  LINE 


AUGUST  1967 
J.  J.  Bussgang 
M.  Leiter 


Prepared  for 

DEPUTY  FOR  SURVEILLANCE  AND  CONTROL  SYSTEMS 
AEROSPACE  INSTRUMENTATION  PROGRAM  OFFICE 

ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
L.  G.  Hanscom  Field,  Bedford,  Massachusetts 


This  document  hos  been  opproved  for  public 
releose  and  sole;  its  distribution  is  un¬ 
limited. 


Project  705B 

Prepared  by 


THE  MITRE  CORPORATION 
Bedford,  Massachusetts 


Contract  AF  19(628)-5  165 


FOREWORD 


The  report  was  prepared  by  the  Range  Communications  Planning  and 
Technology  Subdepartment  of  The  MITRE  Corporation,  Bedford,  Massachu¬ 
setts,  under  Contract  AF  19(628)-5165.  The  work  was  directed  by  the  Devel¬ 
opment  Engineering  Division  under  the  Aerospace  Instrumentation  Program 
Office,  Air  Force  Systems  Command,  Electronics  Systems  Division, 
Laurence  G.  Hanscom  Field,  Bedford,  Massachusetts.  Captain  J.  J.  Centofanti 
served  as  the  Air  Force  Project  Monitor  for  this  program,  identifiable  as 
ESD  (ESSI)  Project  5932,  Range  Digital  Data  Transmission  Improvement. 


REVIEW  AND  APPROVAL 

Publication  of  this  technical  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 


IS  R.  HILL  /Co  lone  1,  USAF 

rector  of  Aerospace  Instrumentation 


Program  Office 


11 


ABSTRACT 


This  paper  considers  differential  phase-shift  transmission  over  a 
band-limited  channel  with  specified  phase  and  amplitude  characteristics. 
The  problem  treated  is  the  prediction  of  feasible  data  rates  and  corre¬ 
sponding  error  rates  over  such  channels. 

An  expression  is  given  for  the  probability  of  error.  This  expression 
takes  into  account  intersymbol  and  interchannel  interference  when  differ¬ 
ential  phase-shift  keying  modulation  is  used.  A  distortion  index  which  is 
a  measure  of  phase  error  is  defined,  and  the  concept  of  phase  opening 
which  parallels  the  concept  of  eye  opening  commonly  employed  in  the 
amplitude-shift  keying  modulation  is  introduced. 

This  analysis  examines  effects  of  various  parameters  on  the  per¬ 
formance  of  a  high  data-rate  transmission  system  over  a  telephone  line. 
Results  indicate  that  for  a  carefully  equalized  line  a  single  tone  system 
should  be  used,  and  that  higher  order  phase  alphabets  are  to  be  preferred 
as  the  bit  rate  increases. 
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SECTION  I 


INTRODUCTION 


In  this  paper,  the  error  performance  of  m-ary  differential  phase- shift 
keying  (PSK)  transmission  of  digital  data  over  toll-grade  voice  channels  is 
considered,  with  special  emphasis  on  high  data  rates.  The  study  takes  into 
account  the  frequency  characteristics  (amplitude  and  delay  distortion  versus 
frequency)  of  the  voice  channel.  It  is  assumed  that  the  channels  are  of  suffi¬ 
ciently  high  grade  that  they  can  support  data  rates  in  excess  of  2400  bits/sec. 
Error  performance  of  the  system  is  analyzed  as  a  function  of  the  number  of 
tones,  the  number  of  possible  phase-shifts  in  each  tone,  the  rate  of  trans¬ 
mission,  the  envelope  shaping  of  the  pulses,  and  the  frequency  separation. 
Several  examples  are  presented. 

In  a  band-limited  channel,  high  data  rates  cause  intersymbol  interference. 
The  effects  of  interference  are  different  for  each  data  sequence.  Thus,  one 
must  compute  the  probability  of  error  for  each  data  sequence  and  find  the 
average  probability  of  error.  The  number  of  sequences  that  must  be  con¬ 
sidered  can  be  so  large  as  to  render  the  time  required  for  the  computation 
of  the  average  probability  of  error  excessive,  even  for  a  high-speed  digital 
computer.  Therefore,  in  addition  to  the  direct  evaluation  of  probability  of 
error,  an  alternate  method  of  evaluating  the  effects  of  the  various  param¬ 
eters  is  examined;  i.e.  ,  the  notion  of  a  distortion  index  introduced  by 
Lucky ^  ^  for  amplitude-shift  keying  is  extended  to  phase-shift  keying. 

This  index  requires  much  less  computation  than  the  probability  of  error 
does,  and  appears  to  be  monotonic  with  error  performance. 
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SECTION  II 


THE  PSK  TRANSMITTER 

One  of  the  m  distinct  characters  is  transmitted  every  T  seconds. 

The  time  interval,  T  ,  and  the  rate  of  character  transmission,  l/T  , 

are  known  as  the  keying  interval  and  the  keying  rate  respectively.  In  a 

phase-shift  keying  system,  the  data  waveform  is  a  segmented  sinewave 

(or  tone)  at  a  frequency  u  /2tt  which  is  an  integer  multiple  of  the  keying 

s 

rate.  The  fundamental  keying  interval  of  duration  T  is  known  also  as  a 

baud  or  a  slot.  The  m  different  characters  correspond  to  the  specified 

differences  in  the  initial  phases  of  the  sinewave  segments  in  successive 

keying  intervals.  Assuming  that  all  the  m  characters  are  equally  probable, 

and  that  the  successive  characters  are  selected  independently,  the  m 

th 

phase-differences  are  chosen  to  be  uniformly  spaced,  i.e. ,  the  f  charac¬ 
ter  corresponds  to  a  phase-difference  27rf/m,  f  =  0,  1,...,  m-1. 

In  order  to  specify  the  transmitted  character,  a  pair  of  successive  sine¬ 
wave  segments  must  be  considered.  To  illustrate  a  pair  of  waveforms,  the 
two  adjacent  keying  intervals  -T/2  <  t  <  T/2  and  T/2  <  t  <  3T/2  are 
arbitrarily  chosen  as  reference  and  datum.  Let  and  be  the  numbers 
corresponding  to  the  transmitted  characters  in  each  of  these  two  intervals 
respectively;  i.e.  ,  the  successive  transmitted  phase-differences  are 
27r  f  ^/m  and  2tt  f  ^/m  .  An  additional  phase-shift  of  tt  /m  radians  is 
purposely  introduced  with  each  keying  interval  in  order  to  minimize  inter¬ 
ference  and  improve  synchronization.  If  this  is  done,  then  the  successive 
transmitted  phases  are  27rf^/m  and  27r(f^  +  f^)/m  +  7r/m  .  In  addition, 
the  transmitted  waveforms  are  assumed  modified  by  an  amplitude  shaping 
function,  k(t),  prior  to  transmission.  Furthermore,  successive  waveforms 
are  allowed  to  overlap  by  6T  seconds  at  each  end,  i.e.  ,  the  sinewave 
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segments  are  (1  +  26)  T  long.  Thus  the  transmitted  waveforms  in  the 
reference  and  datum  intervals  are  represented  by: 

s^(t)  =  k(t)  cos  ^o)gt  +  27rf^/mj 

-  (  \  +  6)  T  *  t  <  +  6)  T  (la) 

s^(t)  =  k(t)  cos  I  a-gt  +  2t r  ^  /m  +  7r/ml 


G  -  6)T  s  1  *  (1  +  6)T  <lb> 

th 

More  generally,  in  the  notation  used  here,  the  r  keying  interval  is 

tVi 

(r  -  1/2)  T  <  t  <  (r  +  1/2)  T  ,  and  the  r  transmitted  waveform  is 
represented  by: 


s  (t) 
r 


=  k(t)  cos 


fw  t 

L  s 


+  2tt 


(*,♦ 


_  1j/m  +  tt  /mj 


(r  -  \  -  6)  T  3  t  -  (r  +  \  *  *)’ 


(2) 


where  2; r  /m  represents  the  phase  transmitted  in  the  interval  r  -  1 

and  f  represents  the  transmitted  character  in  interval  r. 

In  later  sections  of  this  paper  the  Fourier  transform  of  s  (t)  is 
required  and  is  denoted  by  S^(u;)  .  The  computation  of  S^(cr)  can  be 
carried  out  by  inspection  by  treating  each  sinewave  segment  as  a  product 
of  a  rectangular  window  of  duration  (1  +  26)T  and  a  sinewave  of  infinite 
duration. 
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Thus 


S^(o;)  =  F  ^  oj  -  (jJsJ  exp  ^j7r  ^2jC^  +  2L^  ^  +  ij  /m  J 

+  F  ^o>  +  wgj  exp  £  -j7r  (^21^  +  2L^  ^  +  l)  /mj  (3a) 


where 


F(w)  =  ^  sin  [(J  +  fijcoT]  (3b) 


and  K(u;)  is  the  Fourier  transform  of  k(t)  . 
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SECTION  III 


THE  TRANSMISSION  CIRCUIT 

The  transmission  circuit  treated  in  this  report  is  grossly  labeled  as 

the  "telephone  line.  "  The  subscriber  using  a  "telephone  line"  usually  has 

no  clear  description  of  the  physical  circuit.  Most  often  the  actual  circuit 

includes  open  wire  lines,  frequency  multiplex  equipment,  microwave  or 

coaxial  lines,  and  frequently  intermediate  repeaters.  The  properties 

of  such  a  circuit  are  usually  described  by  the  telephone  company  in 

terms  of  bounds  on  the  amplitude  and  delay  distortion  characteristics  in 

the  central  portion  of  the  frequency  region  of  interest.  In  addition,  upper 

bounds  on  the  frequency  translation  error  and  on  noise  hits  are  occasionally 

specified.  Sometimes  a  lower  bound  on  the  expected  signal-to-noise  ratio 

(SNR)  may  also  be  given.  These  specifications  are,  in  fact,  implied  when 

[2] 

the  toll  service  is  described  as,  e.g.  ,  Schedule  4,  Type  A,  B,  or  C. 

The  analyst  concerned  with  the  performance  of  modulation  techniques 

over  "telephone  lines"  is  limited  by  the  imprecise  nature  of  this  description. 

Accurate  predictions  of  performance  are  particularly  difficult  since  echoes 

due  to  even  small  local  irregularities  in  the  frequency  characteristics  are 

[  3] 

known  to  severely  affect  performance.  Since  actual,  precise  charac¬ 
teristics  are  seldom  explicitly  given,  the  analyst  can  generally  draw  only 
certain  broad  conclusions  from  a  few  basic  distortions  of  the  idealized 
characteristics. 

It  is  generally  true  that,  at  conventional  data  rates,  the  contribution 
of  the  thermal  noise  to  the  error  performance  of  a  voice- grade  telephone 
line  is  small,  regardless  of  the  modulation  technique.  At  these  data  rates 
the  performance  of  telephone  circuits  is  limited  by  "impulse  noise"  (due 
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primarily  to  switching  transients)  and  not  by  thermal  noised  "Impulse 
noise"  is  a  term  used  to  characterize  noise  which  has  a  bursty  occurrence 
of  high  amplitude,  short-duration  peaks.  Furthermore,  the  noise  due  to 
crosstalk  from  other  "telephone  lines"  is  usually  insignificant  and  causes 
so  few  errors  that  its  effects  can  be  generally  ignored.  However,  as  the 
data  rate  is  increased,  errors  are  caused  by  intersymbol  and  intertone 
interference,  and  thermal  noise  determines  the  nature  of  the  transition 
from  good  performance  to  poor  performance. 

The  properties  of  the  "telephone  line"  are  represented  by  the  band- 
limited  transfer  function  H(a)  , 


where 

tively, 


H(u>)  = 


\  A  (a; )  exp  j-  j<Mw)  j 


CO 


CO 


CO. 


I  0  elsewhere 

A(a>)  and  <p(u. :)  are  the  amplitude  and  phase  characteristics  respec- 

and  co  and  u>,  are  the  bandwidth  limits  of  angular  frequency, 
a  b 
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SECTION  IV 


THE  RECEIVED  SIGNAL 


The  "telephone  line"  causes  time  spreading  of  the  transmitted  pulse. 

th 

The  output  of  the  telephone  line  in  the  n  keying  interval  due  to  the 

th 

transmitted  sinewave  segment  s  (t)  in  the  r  keying  interval  is  denoted 

r 

by  z  (t)  .  The  Fourier  transform  Z  (u?)  is  readily  found  and  is 
n  ,  r  n  ,  rv  7  J 

simply 


Zn  ,  r(w)  =  [  sr<w)H(w)]  *(JexP  [-  j(n  -  l  )wtJ 

-exp  [-j  ^n  +  J  /jwj  (5) 

where  *  is  the  symbol  for  convolution. 

th 

The  Fourier  transform  of  the  total  received  signal  in  the  n  time 
slot  due  to  an  infinite  string  of  transmitted  sinewave  segments  is  given  by 


Z(w)  « 


I 


Z 

n  , 


r  —  __  OC 

This  summation  for  r  7^  n  represents  the  intersymbol  interference. 


(6) 


In  a  PSK  system  employing  more  than  one  data  tone  on  the  line,  inter¬ 
channel  interference  must  be  included,  in  addition  to  the  intersymbol  interfer¬ 
ence,  i.  e.  ,  the  effect  of  signals  in  all  the  time  intervals  of  other  tone 
frequencies  must  be  considered.  Let  successive  tones  be  numbered  by 
pe(l,  2,  . .  .  M)  and  the  transmitted  tone  designated  by  s.  Let  the  Fourier 
transform  of  the  received  signal  in  the  slot  n  of  the  tone  p  due  to  a 
waveform  in  the  slot  r  of  the  tone  s  be  designated  by  ^(oj  |s). 
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When  p  =  s  this  response  represents  self-interference,  i.  e.  ,  intersymbol 
interference,  and  when  p  ^  s  there  is  intertone  interference. 


If  £  +  L 

r  r 


slot  of  tone  s 
and  imaginary  parts: 


-  1 

then  the  term  Z 


th 

denotes  the  transmitted  phase  in  the  r 
/  \ 

„  Ju  I  S)  can  be  written  as  a  sum  of  real 


ZW  (co  |  s) 
n  ,  r 


Re  jz(P>  (cols)  |  +  jlm  J  Z(p)  (w|s)| 

I  n  ’  r  )  |  n  ,  r  \ 


(7a) 


where 


Re  !z^  (co  1  s)|  =  -I,  (o>)  cos  2n  (  £  +  L 

|  n  ,  r  \  1  '  “ 


r  r  -  1 


/m 


I  (co)  sin  2tt  (  £  +  L  |  /m 

2  V  r  r  -  l) 


Im  \  („|s)[ 

(  n  ,  r  ‘1 


and  where 


cr> 


I’(co)  =  \  do;1  ,  (" K(co'  +  a;)  ±  K(u>’  -  w)  I 

2  J  (jo  —  cu 1  L  J 


CO 


(7b) 


-  I  (a;)  sin  2t  l  £  +  L  ... 

3  \  r  r  -  1  // m 


+  I4(a))  cos  2t  ^  _  1  j/m  (7c) 


sin  ^(co  -  cos  ^  n(a:  -  w')T  +  <K^'')J 


<7d) 
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13  (W)  = 


00  , 


J  co  -  or  L 


K(co  *  +  co )  ±  K  (co '  -  co ) 


CO 


sin 

(c 0  -  co')-y 

sin 

- 

co')T  +  <p{ co') 


(7e) 


The  expressions  for  the  real  and  imaginary  parts  of  Z  ^  (co  |  s)  in  this 

n  ,  r 

form  exhibit  the  two  kinds  of  terms  which  must  be  computed  to  evaluate  the 
performance  of  the  systems:  (1)  the  four  integrals  which  involve  the  fre¬ 
quency  characteristics  of  the  channel  and  the  envelope  shaping  of  the  trans¬ 
mitted  signal,  all  of  which  are  fixed  for  a  given  transmission  system;  and 
(2)  cosine  and  sine  coefficients  which  are  a  function  of  the  transmitted  data 
sequence  alone. 


The  Fourier  transforms  of  the  total  received  data  and  reference  wave¬ 
forms  in  the  n^  and  (n  +  1)  interval  of  the  desired  tone  p  of  an  M 
tone  system  is  given  by 

M  °° 

ZD<“>-  I  I  Zn>,  r<“ls>  (8a> 

s  =  1  r  =  -  00 

M  °° 

z<r  <*»  =  y  y  z<?-i,r<“|s)  <sb> 

s  =  1  r  =  -°° 

The  summation  for  r  is  carried  out  over  both  the  past  and  the  future 
intervals  in  order  to  allow  for  the  case  when  the  specified  idealized 
transfer  function  is  not  physically  realizable. 
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SECTION  V 


THE  PSK  RECEIVER 


The  reason  for  assuming  the  particular  model  of  the  receiver  that  is 
treated  in  this  paper  is  that  this  model  describes  the  principles  embodied 
in  many  current  equipments.  The  receiver  is  designed  to  make,  on  the 
average,  the  fewest  errors  when  the  transmitted  signal  is  perturbed  only 

by  additive  white  Gaussian  noise.  For  each  l£  {l,  2 .  m}  the  receiver 

computes  Pr  {  r^  (t),  r^(t)[i}  ,  i.e.  ,  the  joint  likelihood  of  r^(t)  and 
r^(t)  conditional  on  character  i  having  been  transmitted.  It  is  then 
decided  that  the  transmitted  character  was  that  i  for  which  Pr  (  r^t)  , 
r  (t)  [i}  is  maximum.  Note  that  the  joint  likelihood  is  a  monotonic  function 
X  given  by  Reference  [  5  ]  : 


T/2 


3T/2 


\  dtl  \  dt2rD(tl)rR(t2)coS 

«n(t  -  O  +  2?ri/m 

T/2  T/2 

p  V  1  z  / 

“ 

(9) 


where  is  an  estimate  at  the  receiver  of  the  transmitted  carrier  tone 
P 

frequency  co  . 

The  physical  embodiment  of  the  operations  indicated  by  the  test 

statistic,  X  ,  can  be  obtained  by  three  equivalent  techniques:  (1)  matched 
i 

filtering  and  envelope  detection,  (2)  matched  filtering  and  product  detec- 

[5] 

tion,  and  (3)  multiplication  and  low-pass  filtering . 
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SECTION  VI 


ERROR  PERFORMANCE 

INTRODUCTION  TO  DECISION  PROBABILITIES 

r  g  1 

In  a  previous  paper,  *  the  concept  of  the  "decision  probability”  was 
introduced.  The  decision  probability  is  the  probability  of  choosing  one  out 
of  just  two  (and  not  m)  alternatives.  For  a  PSK  system  this  probability 
is  the  conditional  probability  of  choosing  phase-difference  27rk/m  over 
phase-difference  27ri/m  ,  given  that  an  arbitrary  phase-difference  27rf/m 
was  transmitted  on  tone  p.  For  an  m-phase  system,  this  probability  will 
be  denoted  by  Pr{k  >  i|f  ;  m}  ,  where  the  inequality  sign  indicates  that 
phase-difference  k  is  chosen  over  i.  Decision  probabilities  are  rela¬ 
tively  easy  to  compute  because  they  are  probabilities  of  the  received 
phase-difference  falling  in  half-phases  of  a  phase  diagram.  They  are  of 
interest  because  they  can  be  used  to  determine  the  probability  of  a  charac¬ 
ter  error  which  involves  the  choice  of  one  out  of  all  m  alternatives. 
Moreover,  they  give  the  probability  of  a  binary  digit-error.  Thus  a  two- 
hypothesis  problem  rather  than  an  m-hypothesis  problem  is  solved,  and 
yet  the  desired  m-hypothesis  results  are  obtained.  In  particular,  for 

m-even,  the  probability  of  a  character  error  P  (m)  in  an  m-phase  system 
[51  c 

is  tightly  bound  by 

2Pr  {l  >0|0;m}  -  —  P  (2)  <P  (m)  <  2Pr  {  1  >  0 1 0;m}  (10a) 

1  me  c 

m  >  2 

where 

P  (2)  =  Pr{l  >  0 1 0,2}  (10b) 

c 

is  usually  known  exactly.  For  m  _>  4  ,  the  upper  bound  is  a  good  approxi¬ 
mation  to  Pc(m)  .  Denote  the  probability  of  a  binary  digit-error  for  an 
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m-phase  system  by  P  (m)  ;  for  the  usual  Gray  code  this  probability  can 

[  5 1 

be  expressed  exactly  in  terms  of  the  decision  probabilities,  i.e.  , 

P^(m)  =  Pr{  1  >  0 1  0;m)  m  =  2,4  (lla) 

For  m  =  8  the  probability  P^m)  is  tightly  bounded  by  decision  probabili- 


Pr{  1  >  0 ]  0,8}  +  P  {2  >  1 1 0;8}  -  P  (2) 


<  Pb(8) 


2 

<  - 

3 


Pr{  1  >  0 1  0; 8}  +  Pr  {  2  >  1  ]0;8} 


lib) 


These  results  are  applicable  to  any  maximum  likelihood  choice  between 
m-even  equally-probable  phase-differences.  Similar  expressions  can  be 
derived  for  m-odd  and/or  when  the  phase-differences  are  not  equally- 
probable.  Thus,  computation  of  appropriate  decision  probabilities  permits 
character  error  and  bit  error  probabilities  to  be  evaluated. 


COMPUTATION  OF  DECISION  PROBABILITIES 

The  SNR's  occurring  on  telephone  lines  are  sufficiently  large  so  that 
at  low  data  rates  the  effect  of  thermal  noise  is  to  cause  error  rates 

-7 

typically  below  10  .  Thus  at  low  data  rates,  where  the  distortion  due  to 

the  frequency  characteristic  is  small,  the  impulse  noise  which  causes  error 
-5 

rates  around  10  dominates  over  thermal  noise.  However,  at  high  data 
rates  thermal  noise  does  have  an  effect  on  performance.  It  tends  to  cause 
a  less  abrupt  transition  in  performance  from  the  state  of  no  intersymbol 
or  interchannel  interference  to  the  state  in  which  intersymbol  interference 
is  the  dominant  cause  of  errors.  The  probability  of  error  in  a  telephone 
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line  environment  is  computed  by  taking  advantage  of  the  previous  result  for 

[  5 1 

the  probability  of  error  for  a  channel  perturbed  by  white-Gaussian  noise. 

At  very  high  SNR  (i.e. ,  40  db  or  more,  which  is  typical  for  a  telephone 
line)  a  thresholding  effect  exists,  i.  e. ,  either  relatively  error-free  operation 
or  bottoming  due  to  intersymbol  and  intertone  interference.  The  error  rate 
does  not  improve  further  as  SNR  is  increased. 

In  a  previous  paper  ^  Pr  {  k  >  i  Jf  ;m}  was  obtained  as  a  representa¬ 
tion  of  the  decision  probability  for  a  channel  perturbed  by  additive-white- 
Gaussian  noise  (Ng/2  watt/Hz  for  all  frequencies), 


OO 


i  3  r  1  r  ( 

j  -  2 cos  bs(y  )/N0  +  y  sin  bg(y)  /N0( 

(12a) 


where 


(12b) 


and  where 
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[exp(j2ffi/m)  -  exp(j27r  i/m)]  sin7r(i  +  k  -  21  )/m  I 
cos  2;r  (i  -  f)/m  -  cos2ff(k  -  £)/m  \ 


(12c) 


The  angular  frequency  a?  ,  is  the  estimate  at  the  receiver  of  the  carrier 

P 

frequency  of  the  incoming  keyed  data  tone.  It  is  important  to  note  that 
since  a  and  b  are  a  function  of  the  particular  data  sequence,  the  value  of 
the  integral  in  Equation  (12a)  can  be  different  for  each  data  sequence.  To 
compute  decision  probabilities  exactly,  Pr{k  >  i|f;m}  must  be  computed 
for  all  data  sequences.  An  alternate  approach  will  be  considered  in  the  next 
section. 


[7] 

In  his  analysis  of  2-phase  clairvoyant  PSK,  Lewandowski  transformed 
an  integral  of  the  same  form  as  the  rhs  of  Equation  (12a)  into  a  difference  of 
two  Marcum  Q-functions.  Taking  advantage  of  this  result  Pr|k  >  i|f  ;m) 
can  be  expressed  in  terms  of  Q-functions  as  follows: 


Pr  {k  >  i  £  ;m}  =  1  <1  -  Q 


7a(l  +  b/a)  '  /  a(l  -  b/a)' 

4No  ’  V"1^" 


+  Q 


a(l  -  b/a)  /  a(l  +  b/a) 


4  N 


4  N 


(13a) 


where 


-I 


Q(a  ,  (3)  =  1  -  \  x  IQ  (ax)  exp 


-  (x2  +  a2)/ 2 


dx  (13b) 


While  Equation  (13)  appears  neater  than  Equation  (12),  there  appears  to  be  no 
clear  advantage  to  representing  error  performances  in  terms  of  the  Q-function. 


Unfortunately,  current  tables  of  Q-functions  are  not  generally  available  for 

large  arguments  of  Q  and/or  with  the  accuracy  required  for  the  range  of 

interest  here.  For  example,  a  routine  for  computing  Q-function  with  relative 

-5 

accuracy  AQ/Q  of  10  is  not  sufficient  in  many  problems.  This  can  be 
seen  by  noticing  that 

/a(l  +  b/a)'  j  a(l  -  b/a) 

V  4No  '  V  4No 


in  Equations  (13)  is  approximately  unity  when  Pr{  k  >  i  |f  ;m}  is  small,  so 

_5 

that  a  routine  with  relative  accuracy  AQ/Q  of  10  can  be  used  only  for 
Pr{k  >  i|f;m}  greater  than  10  . 


A  convenient  approximation  for  the  decision  probabilities  is  available 

[  5] 

in  terms  of  the  well-tabulated  error  function. 


Pr{k  >  i|f;m} 


erf 


H1 


) 

\ 


(15a) 


where  erf(x)  is  given  by 


erf(x) 


(15b) 


The  approximation  improves  as  b/a  —  0  .  In  addition  when  b/a  has  the 
maximum  value  of  one,  Equation  (12)  reduces  ^  ^  simply  to,  ^  exp  ^-a/N^ 
The  approximation,  Equation  (15),  is  found  to  be  adequate  for  computations 
considered  in  this  study  when  m  _>  4  ,  but  when  m  =  2  the  exact  expres¬ 
sion,  Equation  (12),  is  required. 
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DISTORTION  INDEX 


In  a  band-limited  channel,  interference  from  symbols  adjacent  in  time 
and  frequency  limits  error  performance.  The  shape  of  the  interference,  and 
hence  its  effect,  depends  on  the  particular  character  sequence.  Hence  in 
order  to  compute  the  average  probability  of  a  binary  digit  error  one  must 
separately  compute  the  probability  of  a  digit  error  for  every  possible  sequence 
of  waveforms.  Since  interference  decreases  with  time  and  frequency  separa¬ 
tion,  in  practice  only  up  to  some  K  intervals  that  cause  interference  above 
the  noise  level  must  be  considered.  Although  each  computation  requires 

only  the  numerical  solution  of  four  integrals,  the  number  of  necessary  compu- 

K 

tations  for  an  m-phase  system  is  m  .  When  the  data  rate  is  high  enough 
compared  to  the  bandwidth,  K  can  become  so  large  that  m  becomes 
enormous.  However,  in  some  cases,  it  may  only  be  desired  to  compare  the 
performance  as  a  function  of  the  number  of  phases/tone,  the  number  of  tones, 
and  the  data  rate.  Therefore,  it  is  useful  to  define  a  quality  index,  which 
would  be  monotonic  with  probability  of  error  and  easier  to  compute  than  the 
probability  of  error  itself. 

The  largest  error  in  the  differential  phase  angle  between  the  reference 
and  data  waveform  that  any  data  sequence  can  produce  was  chosen  as  such 
an  index.  Since  good  (in  terms  of  performance)  sequences  for  some  channels 
may  be  poor  sequences  for  other  channels,  the  problem  of  choosing  a  different 
"worst"  sequence  for  each  channel  remains. 

For  an  amplitude-shift  keying  (ASK)  system,  the  worst  sequence  is  the 
one  having  the  maximum  level  in  all  the  other  keying  intervals.  The  ASK 
distortion  index  is  the  eye  opening  of  the  baseband  signal  at  the  center  of  the 
keying  interval  for  this  sequence.  It  has  been  shown  that  this  index  is 
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monotonic  with  probability  of  error. ^  ^  ^  The  value  of  the  index  is  meaningful 
only  for  channels  and  data  rates  for  which  the  worst  sequence  still  has  an 

open  eye. 

As  an  analogy  for  a  single  tone  PSK  system,  the  "worst"  sequence  is 
defined  as  the  one  whose  characters  selected  in  order  give  the  maxi¬ 
mum  increase  in  phase  error.  The  following  procedure  is  used:  The  phase- 
difference  of  a  pair  of  received  reference  and  data  waveforms  is  computed 
for  each  of  the  m  possible  transmitted  phase  differences.  The  two  symbols 
which  produce  the  maximum  phase  error  due  to  interference  are  retained. 
Then  m  more  computations  are  made  by  considering  the  additional  inter¬ 
ference  due  to  each  of  the  m  possible  phases  in  the  adjacent  slot.  Then 
that  signaling  phase  which  gives  the  maximum  phase  error  for  the  three - 
symbol  sequences  is  retained.  The  procedure  is  continued  slot  by  slot  until 
the  added  interference  due  to  extending  sequences  is  negligible.  Specifically, 


%  =  I“d  ■  “r  -  2lS/m\ 

is  sequentially  computed  where 


a 


R 


-1 

tan 


tan 


Re  Z  p  (*p) 

Z<PD  (“p) 

Re  Z  B  (£p) 

Im  Z  R  (%) 


(16a) 


(16b) 


(16c) 


and  where  2n(. /m  is  the  transmitted  phase  difference.  After  each  computa¬ 
tion  another  slot  is  included  in  ^  and  (u>  ^  while  that 

D  V  P/  R  V  P/ 
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character  sequence  which  maximized  the  previous  D  is  retained.  The 
same  procedure  is  generalized  to  multi-tone  systems. 

Thus,  as  symbols  are  added,  the  already  selected  sequence  is  not 
reconsidered.  Therefore,  there  is  no  guarantee  that  the  selected  sequence 
is  "worst”  in  the  absolute  sense.  However,  the  results  generated  by  this 
procedure  have  been  found  to  be  fully  indicative  of  probabilities  of  error. 
Thus  the  computational  simplification  appears  well  worthwhile. 

In  a  N-level  ASK  system,  the  distortion  index  is  related  to  the 

eye  opening  I  by 

A 


'a  ■  1  -  <N  -  "da 


(17) 


In  a  m-phase  PSK  system,  a  parallel  notion  is  that  of  a  phase  opening 
I ,  which  is  defined  in  terms  of  the  phase  distortion  index  D  .  by 

<P  *  <P  ‘ 


I 


0 


7T 

m 


(18) 


The  angle  1^  represents  a  margin  for  noise  after 
taken  into  account.  The  angle  I  can,  of  course, 

<P 

in  degrees  rather  than  radians. 


intersymbol  interference  is 
be  conveniently  expressed 
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SECTION  VII 


EXAMPLE 

In  order  to  illustrate  the  methods  outlined  in  the  preceding  section, 

consider  the  transmission  characteristics  of  a  nominally  3  kHz  circuit 

shown  in  Figure  1.  These  characteristics  are  representative  of  a  carefully 

equalized  telephone  line.  Let  the  amplitude  shape,  k(t),  be  the  commonly 

used  raised-cosine- shape  overlapping  into  adjacent  slots  on  each  side  by 

one  quarter  of  a  keying  interval  (see  Figure  2).  This  particular  shape  was 

[8] 

found  by  Rappaport  to  be  optimum,  i.e.  ,  to  generate  least  intersymbol 
interference  for  a  quaternary  (m  =  4)  differential  phase  system  at  2400 
bits/sec. 


The  probability  of  error  in  the  binary  data  stream  was  first  computed 
for  2-phase,  4-phase,  and  8-phase  transmissions  on  a  single  tone  located 


Figure  1  .  Characteristics  of  a  Carefully  Equalized  Telephone  Line 
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IB-22,789 


RAISED  COSINE 


Figure  2.  Example  of  an  Envelope  Shaping  Function  k(t) 

at  the  center  of  the  band.  These  results  are  shown  in  Figures  3,  4  and  5. 

It  can  be  seen  that  the  circuit  will  support  2-phase  transmission  at  rates 

up  to  approximately  3600  bit/sec,  4-phase  up  to  6000  bit/sec  and  8-phase 

up  to  7200  bit/sec.  Note  that  40  db  is  a  typical  SNR  and  that  for  2-  and 

-5 

4-phase  transmission,  an  error  rate  better  than  10  is  obtainable.  At 
7200  bit/sec  and  40  db  SNR,  the  8-phase  system  has  an  error  rate  of 

_4 

approximately  10  .  Transmission  of  two  tones  was  next  considered  and 

the  effect  of  frequency  spacing  between  tones  on  the  error  rate  was  examined. 
Results  for  4-phase  DPSK  at  5600  bit/sec  are  plotted  in  Figure  6.  It  can  be 
seen  that  at  this  high  data  rate  the  error  rate  is  very  sensitive  to  proper 
frequency  spacing.  The  optimum  spacing  is  approximately  1450  Hz  which  is 
slightly  above  the  character  rate  of  1400/sec.  This  optimum  point  arises 
from  having  the  band-limited  nature  of  the  channel  counteract  the  inter-tone 
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Figure  3.  Probability  of  Error  versus 
SNR  for  a  1-Tone,  2-Phase  System 
with  a  Raised  Cosine  Envelope 


Figure  4.  Probability  of  Error  versus 
SNR  for  a  1-Tone,  4-Phase  System 
with  a  Raised  Cosine  Envelope 
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Figure  5.  Probability  of  Error  versus 
SNR  for  a  1-Tone,  8- Phase  System 
with  a  Raised  Cosine  Envelope 


Figure  6.  Probability  of  Error  versus 
Tone  Separation  for  a  2-Tone, 

4- Phase  System  at  a  Data  Rate 
Equal  to  5600  bits/sec 
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interference.  When  tones  are  closer  than  1320  Hz  the  intertone  interference 
is  large  enough  to  cause  an  irremovable  phase  error,  hence  "bottoming.  " 

Even  with  the  best  spacing,  a  two  tone  system  was  inferior  to  a  single 
tone  system.  The  same  conclusion  can  be  reached  taking  the  conventional 
"orthogonal"  spacing  (character  rate)  and  examining  the  error  rate  as  the 
number  of  tones  is  increased.  Since  the  overall  data  rate  and  SNR  is  kept 
constant,  as  tones  are  added,  the  keying  rate  and  SNR  of  each  tone  decrease. 
It  has  been  concluded  from  Figure  7  that,  at  a  high  data  rate,  DPSK  over 
this  telephone  line  is  best  operated  with  a  single  tone. 


Continuing  with  the  same  example,  the  phase  opening,  1^  ,  is  plotted  in 
Figure  7  against  the  data  rate  for  2-phase,  4- phase  and  8-phase  single  tone 
systems.  The  computed  curves  indicate  that,  regardless  how  high  the  SNR 
is,  there  is  very  little  phase  margin  for  noise  above  3600  bit/sec  for  a 
2-phase  system,  6000  bit/sec  for  a  4-phase  system,  and  7500  bit/sec  for  an 
8-phase  system.  This  confirms  the  probability  of  error  calculations  in 
Figures  3  and  4  which  are  a  great  deal  more  difficult  to  obtain  than  phase 
opening  calculation.  An  interesting  aspect  of  the  curves  in  Figure  8  is  that 
they  suggest  at  what  bit  rates  one  should  switch  to  a  higher  level  phase- 
alphabet.  For  example,  at  rates  below  3400  bit/sec,  a  2-phase  alphabet 
would  be  recommended.  Between  3400  and  5800  bit/sec  a  4- phase  system 
would  offer  better  performance.  Above  5800  bit/sec  an  8-phase  system 
appears  superior.  Of  course,  at  these  high  data  rates,  the  phase  opening 
is  very  small  and  hence  extremely  large  SNR's  would  be  required. 


In  order  to  verify  that  this  cross-over  in  performance  is  also  exhibited 
by  the  probability  of  error  curves,  2-  and  4-phase  systems  were  compared  at 
data  rates  around  3400  bit/sec  and  SNR's  of  20  db.  At  3200  bit/sec  the  error 

—  6  _ 0 

rates  were  2.0  x  10  and  1.0  x  10  for  4-  and  2-phase  transmission, 
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Figure  7.  Probability  of  Error  versus 
Number  of  Data  Tones 


Figure  8.  Phase  Opening  versus  Data 
Rate  for  a  Raised  Cosine  Pulse 
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respectively.  As  3400  bit/sec,  the  4-phase  error  rate  remained  approximately 
the  same  but  the  2-phase  error  rate  rapidly  deteriorated.  This  confirms  a 
cross-over  slightly  above  3200  bit/sec  for  SNR  =  20  db.  A  higher  SNR  would* 
produce  a  cross-over  at  a  higher  data  rate.  The  cross-over  in  error  rate  per¬ 
formance  around  5800  bit/sec  suggested  by  Figure  8  and  4-  and  8-phase  trans¬ 
mission  can  be  confirmed  by  examining  Figures  4  and  5. 

In  order  to  investigate  what  is  the  effect  of  an  unequalized  transmission 
circuit  on  digital  phase-shift  transmission,  a  time  delay  and  amplitude  charac¬ 
teristic  representative  of  a  Schedule  4  Type  B  circuit  was  assumed  (see 
Figure  9).  The  phase  opening  was  computed  for  2-phase,  4-phase,  and  8- 
phase  transmissions  over  this  circuit  at  different  data  rates.  The  results 
shown  in  Figure  10  indicate  that  2-phase  transmission  could  be  supported 


Figure  9.  Characteristics  Representative  of  a  Schedule  4B  Telephone  Line 
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Figure  10.  Phase  Opening  versus  Data 
Rate  for  a  Raised  Cosine  Pulse 


up  to  rates  slightly  above  3000  bit/sec.  The  interesting  effect,  however,  is 
that  the  higher  order  phase  alphabets  cannot  be  expected  to  increase  the  data 
rates  over  such  a  line. 
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SECTION  VIII 


CONCLUSIONS 

The  results  of  this  study  indicate  that  the  analysis  of  digital  phase- 
shift  transmission  over  band-limited  circuits  can  be  greatly  simplified  by 
computing  the  phase  opening  for  the  worst  data  sequence  rather  than  the 
average  probability  of  error. 

Furthermore,  higher  order  alphabets  can  increase  the  available  data 
rate,  but  only  over  a  well-equalized  line.  Of  course,  at  the  higher  data 
rates,  larger  SNR's  are  required  to  maintain  the  same  error  rate. 


31 


REFERENCES 


1.  R.  W.  Lucky,  A  Functional  Analysis  Relating  Variation  and  Intersymbol 
Interference  in  Data  Transmission,  BSTJ,  42,  5,  (Sept.  ,  1963),  2427-2484. 

2.  American  Telephone  and  Telegraph  Co.  ,  Tariff  FCC  #260,  22  May 
1966,  p.  65. 

3.  H.  A.  Wheeler,  The  Interpretation  of  Amplitude  and  Phase  Distortion 
in  Terms  of  Paired  Echoes,  Proc.  IRE,  27,  (June,  1939),  359-385. 

4.  A.  B.  Bodonyi,  Effects  on  Impulse  Noise  on  Digital  Data  Transmission, 

IRE,  CS-9,  (Dec.,  1961),  355-361. 

5.  J.  J.  Buss  gang  and  M.  Leiter,  Analysis  of  a  General  Class  of  Phase- 
Shift  Keying  Communication  Systems,  Cruft  Laboratories,  Technical 
Report  No.  494,  Harvard  University,  Cambridge,  Mass.,  Feb.,  1966. 

6.  J.  J.  Bussgang  and  M.  Leiter,  Error  Rate  Approximations  for 
Differential  Phase-Shift  Keying,  IEEE  Trans.  ,  CS-12,  (Mar.  ,  1964), 

18-27. 

7.  L.  Lewandowski,  Derivation  of  the  Probability  of  Error  in  a  PSK  System 
Employing  a  Phase  Locked  Loop,  Internal  Memorandum,  Syl vania 
Amherst  Laboratories,  Syl  vania  Electronic  Systems- Central,  Buffalo, 

New  York,  7  Feb.  ,  1962. 

8.  M.  A.  Rappaport,  Digital  Computer  Simulation  of  a  Four-phase  Data 

Transmission  System,  BSTJ,  43,  3,  (May,  1964),  927-964. 


32 


Unclassified _ 

_ S«h  urity  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 

(S^inrify  c  la  s  m  fu at  ion  of  title,  body  of  abstract  and  indexing  annotation  must  t>c  entered  wh  rn  the  overall  report  is.  classified) 


i  ORIGINATING  ACTIVITY  ((’ofporAfe 

The  MITRE  Corporation 

Bedford,  Massachusetts 

2a.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

2b.  GROUP 

3  R  1  PORT  TITLE 

Error  Performance  of  Differential  Phase  Shift 

Transmission  over  a  Telephone  Line 

4  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  dates) 

n/a 

s  au  THORtSl  (First  name,  middle  initial,  last  name) 

Bussgang,  Julian  J. ;  Leiter,  Myron 

6  REPOR  T  O  A  TE 

August  1967 

7a.  TOTAL  NO.  OF  PAGES  7 b.  NO.  OF  REFS 

37  8 

8*.  CONTRACT  OR  GRANT  NO 

AF  19(628)-5165 

b  PROJ  EC  T  NO 

705B 

C. 

d. 

9a.  ORIGINATOR'S  REPORT  NUMBER!?) 

ESD-TR -67-393 

9b.  OTHER  REPORT  NO(S>  (Any  other  numbers  that  may  be  assigned 
this  report) 

MTP-66 

10  DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 

11  SUPPLEMENTARY  NOTES 

12.  SPONSORING  MILITARY  ACTIVITY  DgpULy  for  SUFVeil- 

lance  and  Control  Systems,  Aerospace  Instru¬ 
mentation  Office;  Electronic  Systems  Division, 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

13  ABSTRACT 


This  paper  considers  differential  phase-shift  transmission  over  a  band-limited  channel 
with  specified  phase  and  amplitude  characteristics.  The  problem  treated  is  the  prediction  of 
feasible  data  rates  and  corresponding  error  rates  over  such  channels. 

An  expression  is  given  for  the  probability  of  error.  This  expression  takes  into  account 
intersymbol  and  interchannel  interference  when  differential  phase-shift  keying  modulation  is 
used.  A  distortion  index  which  is  a  measure  of  phase  error  is  defined,  and  the  concept  of 
phase  opening  which  parallels  the  concept  of  eye  opening  commonly  employed  in  the  amplitude- 
shift  keying  modulation  is  introduced. 

This  analysis  examines  effects  of  various  parameters  on  the  performance  of  a  high  data- 
rate  transmission  system  over  a  telephone  line.  Results  indicate  that  for  a  carefully  equalized 
line  a  single  tone  system  should  be  used,  and  that  higher  order  phase  alphabets  are  to  be 
preferred  as  the  bit  increases. 
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